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Abstract 

We calculate the A-dependence of charm and beauty production cross sections 
on nuclear targets taking into account the difference of quark and gluon distribu- 
tions in free nucleons and in nucleus. At comparatively low energies, if a ~ A a , 
the value of a is slightly higher than unity. With the growth of the initial energy 
the value of a decreases and becomes smaller than unity. We also calculate the 
values of a for different Feynman-x of the produced QQ pair and obtain that they 
decrease significantly in the beam fragmentation region. 
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1 Introduction 



The processes of heavy flavour production on nucleon and nucleus targets at high ener- 
gies are very interesting from both the theoretical and practical points of view. These 
processes at high energies are usually considered in the framework of perturbative QCD 
in the leading and the next-to-leading order a s expansion. In the case of a nuclear target 
it is usually assumed that the cross section of heavy flavour production, a(QQ), should 
be proportional to A a , with a — 1, that is, in agreement with the experimental result of 
[1], where a value a = 1.02 ± 0.03 ± 0.02 was obtained for the case of neutral D-meson 
production at y/s = 39 GeV. However it is experimentally well-known [2-8] that quark 
distributions in nucleus are slightly different from the same distributions in free nucle- 
ons. So the value of a can differ from unity and it seems interesting to estimate it. We 
will present in the Section 2 the results obtained using the nucleus to nucleon structure 
function ratios taken from Ref. [9]. 



2 QCD predictions of charm and beauty production 
A-dependence 

The standard QCD expression for heavy quark production cross section in a hadron 1 - 
hadron 2 collision has the form 



a 



12— QQ 



X a Ga/l(x a , p?) XbGb/2(Xb, /i 2 ) 



a ab ^(s,m Q ,p 2 ), (I) 



where x a o = Atrq/s and Xbo = 4mq/ (sx a ). Here G a /i(x a , p 2 ) and Gb/2(xb, p 2 ) are the 
structure functions of partons a and b inside hadrons 1 and 2 respectively, and 
fr ab -*QQ(s,mQ, fi 2 ) is the cross section of the subprocess ab — > QQ as given by standard 
QCD. The latter depends on the parton center-of-mass energy s = (p a + pb) 2 = x a XbS, 
the mass of the produced heavy quark tuq and the QCD scale /x 2 . Eq. ([!]) should account 
for all possible subprocesses ab — > QQ. 

The parton cross section cr ab ^^^(s,mQ, fi 2 ) can be written in the form |10 



a 2 (/i 2 ^ 



a ab ^(s,m Q ,n 2 ) 



m 



2 fab{p^ 2 /m 2 Q ), 



with 
and 



P = 4m 2 Q /s 

f ab (p, p 2 /m 2 Q ) = f%\p) + Ana s (p 2 )[fH\p) + ff b \p) In (/i 2 /m 2 )] 



(2) 

(3) 
(4) 



2 



The functions f^(p), fil\p) an d fal\p) can be found in [|i~0f 



For the numerical calculations we wrote the nuclear structure function G b /A_{x b , P 2 ) in 
the form 

G b/A (x b , /j 2 ) = A- G b/N (x b , /j 2 ) ■ R b (x, /i 2 ), (5) 

similarly to Ref. [|Il|] and we take the values of R^x, /i 2 ) for gluons, valence and sea 
quarks from Ref. [9]Q. They are presented in Fig. 1. The values of R^(x,fi 2 ) in [9] 
are presented for x > 10 -3 that is not small enough at high energies. So at x < 10~ 3 
we used two variants of the R b (x, /i 2 ) behaviour for gluon distributions: The first is the 
extrapolation as x^ with f3 = 0.096 and 0.040 for the charm and beauty production, 
respectively (solid lines in Fig. 1). Such behaviour is in qualitative agreement with the 
results of fL3| . The second is the constant frozen at x = 10~ 3 (dash-dotted lines). Three 
different sets of parton distributions were used, namely MRS-1 ||14|| , MT S-DIS [15] and 
GRV HO QT§ that can be found in CERN PDFLIB 0, but they give practically the same 
results for the a behaviour (see Table 1). We have used in the case of charm production 
the values m c = 1.5 GeV and /i 2 = 4 GeV 2 and in the case of beauty production m b = 
5 GeV and /z 2 = m 2 . 

The obtained results for a determined from the ratios of heavy quark production cross 
section on a gold target and on the proton are presented in the Table 1 for three sets 
of structure functions and frozen gluon distribution ratios at x = 10 -3 and in Fig. 2 
for the GRV HO set and two variants of gluon distribution ratios in the small x region. 
Here ^/snn is the cm. energy for the interaction of the incident proton with one target 
nucleon. One can see that at fixed target energies the values of a are slightly higher than 
unity, which is not in contradiction with the result of Ref. [1] . However a decreases with 
increasing energy and this effect is larger in the case of charm production than in the 
case of beauty. One can see also that the difference between two variants for i?^(x,/x 2 ) 
at x < 10 -3 becomes important only at the highest energies. 

We also calculate the values of a for different Feynman-x (xf) regions using xf = 
x a — x b at energies y/s^N = 39 GeV and 1800 GeV. The results are presented in Fig. 
3. At negative and moderate xp (in the nucleus fragmentation region) the values of a 
are slightly higher than unity. However in the case of charm production in the beam 
fragmentation region (positive xf) the values of a become essentially smaller than unity. 
For beauty production the last effect is expected only at very high energies. 



^et us note that the theoretical estimations in an earlier paper |12| give more or less a similar 
behaviour of i?^(x,/i 2 ) for light sea quarks but an enhancement of gluon distributions in the nuclei at 
small x. 
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3 Conclusions 



We calculate the A-dependence of charm and beauty production using standard QCD 
formulas and accounting in the difference of parton distributions for free and bound 
nucleons. If one parametrize the heavy flavour production cross section as o ~ A a , the 
value of a is slightly different from unity at the available energies. At comparatively low 
energies the obtained values of a are a little larger than unity This should be connected 
with some nucleon-nucleon correlations which change the large- x parton distributions. 

At higher energies the values of a decrease and become smaller than unity. At ^/snn 
= 1800 GeV we expect the value of a ~ 0.95. The decrease of the ratios R^x,^ 2 ) that 
results in a decrease of a can be connected with the effects of parton density saturation 



18| which occurs in heavy nuclei at x values higher than in the proton. 



If we consider two small and different values of x a and Xj, in Eq. (1), it is clear that 
the contribution to the inclusive cross section from the region x a < xj, should be larger 
than the mirrow contribution (x a > Xb) because the value of the ratio R^(x, /i 2 ) in the 
first case is larger. It means that heavy quark pairs will be produced preferably in the 
nucleus fragmentation hemisphere, i.e., asymetrically, that is quite usual and confirmed 
experimentally in the case of light quark hadron production. From Fig. 3 it is clear that 
charm production on nucleus at LHC energies should give information on the nuclear 
shadowing of the structure functions at small x. 

Let us note also that at HERA-B energy (y/s^N — 39 GeV) with GRV(HO) parton 
distributions we obtain a(pp — ► bb) = 6.8 nb that is in agreement with the experimental 
value 5.7± 1.5 ± 1.3 nb of Ref. ||19|| . However the experimental point should be decreased 
by 5% because it was obtained as an extrapolation from the gold target data assuming 
an A 1 dependence whereas our calculations give the value a = 1.01. 

We are grateful to M.A.Braun for useful discussions and both to K.J.Eskola and to 
I.Sarcevic for sending us their numerical results. We thank the Direction General de 
Polftica Cientifica and the CICYT of Spain for financial support. N.A. and C.A.S. also 
thank the Xunta de Galicia for financial support. The paper was also supported in part 
by INTAS grant 93-0079. 
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Table 1 



Perturbative QCD predictions for the a values of charm and beauty production in high 
energy pA interactions with gluon distribution ratios frozen at x — 1CT 3 . 



y/s, GeV 


MSR-1 


MT S-DIS 


GRV(HO) 


cc 


bb 


cc 


bb 


cc 


bb 


27 


1.01 


1.01 


1.01 


1.00 


1.01 


1.01 


39 


1.01 


1.01 


1.01 


1.01 


1.01 


1.01 


62 


1.00 


1.01 


1.00 


1.01 


1.00 


1.01 


120 


0.99 


1.01 


0.99 


1.01 


0.99 


1.01 


200 


0.99 


1.01 


0.99 


1.01 


0.98 


1.00 


1800 


0.96 


0.99 


0.96 


0.99 


0.95 


0.99 


14000 


0.94 


0.98 


0.94 


0.98 


0.93 


0.97 
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Figure captions 



Fig. 1 : The functions Rq(x,^ 2 ), Ry(x, /i 2 ) and Rg(x,fi 2 ) which determine the ratios of 
the distributions for protons in the nucleus versus free protons for gluons (solid 
and dash-dotted curves, see text), valence quarks (dashed curves) and sea quarks 
(dotted curves) respectively for fj, 2 = 4 GeV 2 (a) and /i 2 = m 2 (b). 

Fig. 2 : The energy dependence of a in the cases of charm and beauty production for 
GRV HO structure functions and using extrapolated (solid curves) and frozen 
at x = 1CT 3 (dashed curves) ratios of gluon distributions. 

Fig. 3 : The Feynman-x dependence of a values in the cases of charm and beauty produc- 
tion at ^snn = 39 GeV (a) and 1800 GeV (b) for GRV HO structure functions 
and using extrapolated (solid curves) and frozen at x — 10~ 3 (dashed curves) 
ratios of gluon distributions. 
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Fig . 1 a 
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